Naval  Research  Laboratory 

Washington,  OC  20375-5000 


NRL  Report  9280 


< 

i 

D 

< 


A  Framework  for  the  Detection  of  a  Target 
of  Unknown  Velocity 

M.  J.  Steiner 

Target  Characteristics  Branch 
Radar  Division 


September  27,  1990 


90  10  23  152 

Approved  for  public  release;  distribution  unlimited 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704  0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  i  nour  D*r  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this 
collection  of  information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services.  Directorate  for  information  Operations  and  Reports,  12 IS  Jefferson 
Davis  Highway,  Suite  1204.  Arlington,  va  22202-4302.  and  to  the  Office  of  Management  and  Budget.  Paperwork  Reduction  Project  {0704-0188).  Washington,  DC  20503 


1.  AGENCY  USE  ONLY  (Leave  blank) 


2.  REPORT  DATE 

September  27,  1990 


4.  TITLE  AND  SUBTITLE 

A  Framework  for  the  Detection  of  a  Target 
of  Unknown  Velocity 


6.  AUTHOR(S) 

Steiner,  Michael 


3.  REPORT  TYPE  AND  DATES  COVERED 

Interim 


S.  FUNDING  NUMBERS 


PE  -  PE61153N 
TA  -  TARR02 1-05-13 
WU  -  DN480-006 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADORESS(ES) 

Naval  Research  Laboratory 
Washingon,  DC  20375-5000 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

NRL  Report  9280 


9.  SPONSORING /MONITORING  AGENCY  NAME{S)  AND  AODRESS(ES) 

Office  of  Naval  Research 
Arlington,  VA  22217 


10.  SPONSORING  /MONITORING 
AGENCY  REPORT  NUMBER 


12a.  DISTRIBUTION /AVAILABILITY  STATEMENT 


12b.  DISTRIBUTION  CODE 


Approved  for  public  release;  distribution  unlimited. 


ABSTRACT  (Maximum  200  words) 


The  detection  of  a  target  in  correlated  Gaussian  clutter,  thermal  noise,  and  extraneous  interference  is  con¬ 
sidered.  The  amplitude,  phase,  and  Doppler  frequency  of  the  signal  are  not  known  a  priori.  The  optimal  detec¬ 
tor,  in  the  Neyman-Pearson  sense  is  difficult  to  implement,  and  suboptimal  methods  are  used.  A  general  cri¬ 
terion  that  measures  the  performance  of  suboptimal  detectors  relative  to  an  optimal  test  is  presented.  This  cri¬ 
terion  is  encompassed  in  a  design  procedure  used  to  design  Doppler  filters.  The  procedure  allows  many  design 
considerations  to  be  taken  into  account,  which  results  in  a  design  that  attempts  to  minimize  the  number  of  filters 
required.  For  low  dimensionality  the  procedure  is  of  single  filter  designs;  for  higher  dimensionality  multiple 
filters  are  designed.  The  performance  of  these  is  compared  with  the  results  obtained  by  Emerson.  Andrews, 
and  the  generalized  likelihood  ratio  test.  The  clutter  to  thermal  noise  ratio  and  spectral  width  of  the  clutter  are 
assumed  known.  When  these  parameters  are  not  known  exactly,  the  robustness  of  these  and  other  filter  systems 
is  examined.  A  robustness  theorem  is  presented,  which  relates  the  performance  of  a  fixed  filter  to  the  clutter  to 
thermal  noise  ratio  as  a  function  of  the  Doppler  frequency.  Future  research  is  discussed.  s  t 


14.  SUBJECT  TERMS 

Detection 
Space-based  radar 


17.  SECURITY  CLASSIFICATION 
OF  REPORT 

UNCLASSIFIED 


MSN  7540-0 1  -280-5500 


Filter  design 


18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

UNCLASSIFIED 


SECURITY  CLASSIFICATION 
OF  ABSTRACT 

UNCLASSIFIED 


IS.  NUMBER  OF  PAGES 

38 


16  PRICE  CODE 


20.  LIMITATION  OF  ABSTRACT 


Standard  Form  298  (Rev  2  89) 

P'PSCNWKl  by  ANSI  Std  239  ’8 
298  102 


CONTENTS 


1.  INTRODUCTION  .  1 

2.  DEFINITIONS  .  2 

3.  MINIMAX  LOG  ENERGY  CRITERIA  (MMLE)  .  4 

4.  MMLE  APPLIED  TO  THE  RADAR  PROBLEM  .  6 

5.  DESIGN  FOR  LOW-DIMENSION  N  .  7 

6.  DESIGN  FOR  HIGH-DIMENSION  N .  13 

7.  ROBUSTNESS  .  17 

8.  CONCLUSIONS  . 23 

REFERENCES  . 24 

APPENDIX  A  —  Improvement  Factor  . 25 

APPENDIX  B  —  Proof  of  Robustness  Theorem  . 29 

APPENDIX  C  —  Computation  of  MMLE  Filters  .  33 


teeesslon  For  y\ 

]  NTIS  GRAM 
DTIC  TAB 
i  Unannounced 
]  Justificatio 

;] 

Wf 

□ 

□ 

n 

i 

'  Ry 

Distribution/ 

Aval 

Dlst 

Labilit 

Avail 

Spec 

y  Codas 
ind/or 

Lai 

iii 


A  FRAMEWORK  FOR  THE  DETECTION  OF  A  TARGET 
OF  UNKNOWN  VELOCITY 


1.  INTRODUCTION 

A  basic  problem  in  detection  theory  is  that  of  detecting  a  target  of  unknown  velocity.  In  each 
resolution  cell  of  a  pulsed  radar,  a  detector  is  applied  to  determine  whether  interference  only  is 
present  ( H0 )  or  a  (target)  signal  is  also  present  (H ,).  The  interference  is  composed  of  the  sum  of 
thermal  noise  and  clutter.  Generally  three  signal  parameters  are  not  known  a  priori:  amplitude,  phase, 
and  Doppler  frequency.  When  the  distribution  of  these  variables  is  known  a  priori,  the  resulting 
hypothesis  test  can  be  solved  through  the  Neyman-Pearson  test.  In  this  case,  a  likelihood  ratio  is  aver¬ 
aged  over  the  random  parameters  [1,2].  Although  the  resulting  test  satisfies  the  theoretical  require¬ 
ments  of  optimality,  its  complexity  renders  it  too  difficult  to  implement.  Suboptimal  methods  are 
employed  to  approximate  the  performance  of  the  average-likelihood  ratio. 

The  generalized  likelihood  ratio  test  (GLRT)  combines  maximum  likelihood  estimates  of  the 
unknown  parameters  with  the  likelihood  ratio  [3].  This  is  discussed  in  Section  2,  Eq.  (6)  for  the  radar 
problem.  Essentially  Eq.  (6)  requires  finding  the  maximum  over  frequency  of  the  magnitudes  of  the 
outputs  of  linear  filters.  Such  a  test  is  also  difficult  to  implement  because  the  maximization  is  over  all 
possible  values  of  frequency.  Suppose  the  frequencies  over  which  the  maximization  takes  place  are 
partitioned  into  a  finite  discrete  set  {/i ,  /2,  •  •  •  ,  //>},  the  resulting  test  (finite  constrained  GLRT  or 
CGLRT)  now  requires  evaluation  of  only  P  linear  filters.  If  P  is  low,  the  implementation  requires 
substantially  less  hardware  compared  to  the  unconstrained  test  or  GLRT.  However,  such  a  test  may 
fail  to  perform  close  to  what  can  optimally  be  achieved.  Under  such  a  constraint  on  the  number  of 
filters,  it  is  necessary  to  find  alternative  means  for  designing  filters  that  perform  close  to  what  can 
optimally  be  achieved.  When  the  number  of  retuns  N  is  small,  a  single  filter  may  suffice  to  approxi¬ 
mate  the  performance  of  the  cumbersome  average-likelihood  test.  These  single  filters  are  commonly 
referred  to  as  moving  target  indicators  (MTIs).  Heuristically,  the  clutter  is  notched  around  zero 
Doppler,  while  moving  targets  and  thermal  noise  are  passed.  For  larger  N,  single  filters  are  not  as 
effective  as  multiple  filters  that  can  also  integrate  the  target  out  of  the  thermal  noise. 

The  design  of  single  nonrecursive  filters  for  low  N  is  presently  based  on  a  technique  first 
developed  by  Emerson  and  later  expanded  by  Capon  [4,5].  This  procedure  minimizes  the  interfer¬ 
ence  output  power  or  equivalently,  maximizes  the  average  signal  to  interference  (SIR)  ratio  at  the  out¬ 
put  of  the  filter  when  the  target  Doppler  shift  is  uniformly  distributed  over  all  values  [6],  The 
interference  consists  of  the  sum  of  thermal  noise  and  clutter  and  typically  is  dominated  by  the  clutter. 
The  eigenvector  corresponding  to  the  minimum  eigenvalue  of  the  covariance  matrix  of  the  interfer¬ 
ence  is  the  weight  to  be  used  in  the  filter.  The  interference  output  energy  of  such  a  (normalized) 
filter  is  the  minimum  eigenvalue  of  the  covariance  matrix  of  the  interference.  Andrews  [7]  has 
solved  the  design  for  multiple  filters  based  on  this  criterion. 

Although  Emerson’s  criterion  works  well  for  many  typical-interference  power  spectra,  it  fails 
when  the  power  spectra  are  modified  in  other  applications.  The  problems  take  the  form  of  holes  in  the 
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Doppler  coverage  whereby  certain  target  speeds  cannot  be  detected.  This  is  illustrated  in  Appendix 
A.  Additionally,  Emerson’s  procedure  is  unable  to  provide  a  desired  target  velocity  response  weight¬ 
ing.  This  latter  problem  has  been  addressed  by  modifying  Emerson’s  filter  by  quadratic  program¬ 
ming  [8], 

The  limitations  seen  in  Appendix  A  are  inherent  to  the  form  of  the  criteria  as  discussed  in  Sec¬ 
tion  2.  The  purpose  of  this  report  is  to  present  a  criterion  that  does  not  have  these  limitations  and 
produces  filters  with  desirable  characteristics.  This  criterion  is  presented  in  Section  3.  In  Section  4 
the  criterion  is  applied  to  the  radar  problem,  and  an  equation  is  derived  for  the  design  of  filters.  This 
result  is  used  in  Section  5  to  design  single  filters  for  small  N.  Comparisons  are  made  with  Emerson’s 
method.  In  Section  6,  multiple  filters  are  constructed  and  the  results  are  compared  with  Andrews’ 
method  and  the  CGLRT.  In  Section  7,  the  robustness  of  filters  is  examined  relative  to  the  clutter  to 
thermal  noise  ratio  C/N0  and  the  normalized  spectral  width  of  the  clutter.  The  results  are  summar¬ 
ized  in  Section  8. 


2.  DEFINITIONS 


The  classical  hypothesis  testing  problem  encountered  when  trying  to  discriminate  a  moving  tar¬ 
get  from  interference  is  represented  as 


Hq  :  Y  =  X 


1 

ej8 


(1) 


Hx  :  Y  =  X  +  seJ*Ss  >0  S  = 


ej(N  -  iw 
V-  J 

where  N  denotes  the  dimension.  The  signal  vector  is  S.  The  vector  X  consists  of  the  sum  of  thermal 
noise  and  clutter  and  is  assumed  to  have  zero  mean  and  a  Gaussian  distribution.  A  covariance  matrix 
R  associated  with  X  is  defined  as  R}  k  =  E(xjxk),  where  *  denotes  complex  conjugate.  The  diagonal 
elements  of  R  are  fixed  at  a2,  which  is  the  sum  of  the  thermal  and  clutter  energy  per  component. 
The  unknown  constant  phase  <t>  of  the  received  signal  is  assumed  to  have  a  uniform  distribution  in  the 
interval  [0,  27r],  The  Doppler  frequency  6  is  also  assumed  random  with  a  uniform  distribution  in  the 
interval  [0,  27t].  The  signal  amplitude  s  is  assumed  to  be  greater  than  zero;  the  prior  distribution  of  s 
is  unknown.  The  probability  of  Y  given  s,<t>,6  assuming  the  thermal  noise  is  Gaussian  distributed  is 


P(Y\s,4>,9) 


1  -<Y-.«>',>S )’R  '(Y-  sp'*S)* 

(7 r)N  I  R  | 


(2) 


where  T  denotes  vector  transpose.  The  likehood  ratio  is. 


L(Y  |  s,<j),8) 


P,(Y|j.<M) 

/VY) 


(3) 


The  optimal  test  as  discussed  in  Section  1  is  the  Neyman-Pearson  or  average-likelihood  ratio  test 
[9)  that  compares  the  average-likelihood  with  a  threshold.  The  average  likelihood  can  be  written  as 
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L{ Y  |  i) 


f  *  1°°  L(Y|  s,4>,d)P(d>)  P($)d<t>dd 

*  _i‘sr*  'S70(2s  I  Y^-'S*  \)~d6. 

U  7  nr 


(4) 


This  detector  is  difficult  to  implement.  The  test  can  be  approximated  and  considerably 
simplified  by  using  the  GLRT  that  compares 

maxL(Y  |  s,<t>,6)  (5) 

i.  <t>.6 

to  a  threshold.  This  can  be  simplified  for  our  problem  to 


1 


6(Y)  = 


<  0  or  1  if  max0  K(d )  |  Y TR  'S*  | 
0 


> 

=  T 
< 


(6) 


where  K(d)  =  l/^ST  R~l S*.  This  test  has  a  number  of  desirable  properties  when  the  frequency  0 
is  fixed,  (1)  it  can  be  accomplished  by  taking  the  magnitude  of  a  linear  filter,  (2)  it  is  equivalent  to 
the  average  likelihood  test  in  Eq.  (4)  hence  is  optimal  for  fixed  frequency,  and  (3)  it  possesses  the 
desirable  property  of  being  uniformly  most  powerful  (UMP)  with  respect  to  the  unknown  amplitudes. 
Although  methods  exist  for  finding  the  maximum  of  the  above  equation  over  0,  they  are  difficult  to 
implement.  Estimation  methods  for  finding  0  for  a  similar  problem  are  given  in  Ref.  10.  However, 
if  |  \tR  -iS*  |  is  implemented  for  each  0,  a  matched  filter  needs  to  be  built  or  evaluated  digitally  for 
each  0.  This  requires  an  infinite  number  of  filters.  Therefore  the  test  is  typically  approximated  by  a 
finite  number  of  filters.  In  this  case  each  filter  must  cover  an  interval  of  frequencies  0.  Degradation, 
relative  to  the  infinite  bank  of  filters,  Eq.  (6),  will  occur  at  frequencies  where  the  finite  filters  become 
mismatched.  An  optimality  or  design  criterion  to  measure  this  degradation  can  be  defined,  and  new 
coefficients  can  be  determined.  The  amount  of  the  deterioration  with  frequency  can  vary  greatly, 
depending  on  the  filters  used,  the  covariance  matrix,  and  on  the  length  N.  The  best  filter  design  when 
the  number  of  filters  is  limited  along  with  the  associated  design  criteria  is  paramount  to  performance 
the  cost. 


An  optimal  criterion  that  results  from  hypothesis  testing  is  the  maximum  of  the  average  proba¬ 
bility  of  detection  (Pd)  over  the  frequency  range.  For  independent  and  identically  distributed  (HD) 
Gaussian  interference  with  a  Raleigh  distribution  on  the  amplitudes,  Brennan  et  al.  [1]  show  that  the 
average  likelihood  ratio  test  (which  maximizes  the  average  probability  of  detection)  and  the  general¬ 
ized  likelihood  ratio  test  perform  in  nearly  the  same  manner.  Little  has  been  done  to  design  single 
filters  when  the  criterion  of  maximizing  the  average  probability  of  detection  is  desired  when  corre¬ 
lated  interference  is  present.  In  fact  this  is  not  a  good  criterion.  Since  the  problem  of  correlated 
Gaussian  interference,  Pd  may  range  over  several  orders  of  magnitude  depending  on  the  frequency  of 
the  target.  Hence  an  average  of  these  values  over  the  frequency  range  is  strongly  dominated  by  the 
largest  values.  The  frequencies  where  the  largest  values  occur  may  be  over  a  narrow  interval  that  may 
not  even  be  considered  important.  The  second  criterion  that  predominates  in  the  literature  is  the  one 
of  Emerson,  which  consists  of  minimizing  the  average  interference  power  at  the  output  of  the  filter. 
The  problems  associated  with  this  criterion  are  discussed  in  Appendix  A. 
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An  optimal  criterion  is  needed  that  can  be  used  to  design  either  a  single  filter  or  a  bank  of  filters 
that  do  not  suffer  the  limitations  of  the  previous  methods.  The  properties  that  we  consider  necessary 
for  the  criterion  to  satisfy  are  the  following: 


•  The  criterion  should  accurately  reflect  the  degradation  of  the  desired  filters  relative  to  what 
can  be  achieved  in  an  optimal  sense; 

•  The  criterion  should  be  broad  enough  to  allow  parameters  of  the  design  (such  as  number  of 
filters,  positioning  of  the  filters,  etc.)  to  be  varied; 

•  The  criterion  should  be  capable  of  being  weighted  to  take  a  particular  desired  profile  or  cost 
function  into  consideration. 

For  the  criterion  to  reflect  the  degradation  of  the  bank  of  filters,  a  suitable  distance  measure 
must  be  defined.  An  appealing  indicator  at  a  particular  frequency  0  is  Pd  for  a  given  probability  of 
false  alarm  (P/a).  Obviously  we  would  like  Pd  to  be  as  close  to  what  can  be  optimally  achieved  at  a 
given  a  Pfa.  To  determine  what  values  are  actually  close,  an  energy  function  can  be  introduced.  This 
function  should  tell  us  how  much  extra  energy  is  needed  for  a  suboptimal  processor  to  perform  as 
well  as  an  optimal  processer  as  a  function  of  frequency.  Recently,  Compemolle  [11]  has  brought  to 
attention  a  distance  measure  for  use  in  speech  recognition  for  measuring  the  difference  between  power 
spectrums.  This  measure  is  based  upon  the  difference  of  the  two  power  spectra  in  the  log  domain. 
Similarly,  the  radar  engineer  is  interested  in  the  signal  energy  and  noise  energy  expressed  in  deci¬ 
bels.  The  reason  for  preferring  decibels  over  absolute  energy  values  relates  to  the  radar  range  equa¬ 
tion  and  the  exponential  nature  of  the  probability  of  error  equations.  A  more  relevant  distance  meas¬ 
ure  is  based  on  the  log  energy.  In  effect,  it  is  desirable  to  know  as  a  function  of  0  the  necessary  sig¬ 
nal  energy  in  dB  required  of  the  suboptimal  processor  to  perform  as  well  as  the  optimal  processor. 
A  criterion  based  on  this  measure  is  introduced  and  used  in  the  remainder  of  the  report. 

3.  MINIMAX  LOG  ENERGY  CRITERIA  (MMLE) 

A  general  optimal  criterion  that  satisfies  the  properties  we  required  of  an  optimal  criterion  is 
now  presented  in  Section  2.  The  input  SIR  s2 /a2  is  denoted  as  7. 

Suppose,  for  each  frequency  0  €  [0,  2tc]  the  Neyman-Pearson  design  criteria  is  used  to  find  the 
most  powerful  tests  5?(Y),  where  7  denotes  the  input  SIR,  under  the  constraint  Pfa(5$(Y))  =  a. 
Then  any  rule  S?(Y)  will,  by  the  Neyman-Pearson  lemma,  satisfy 

WOO !  9  =  6)  <  Pd(6e 7(Y)  1  9  =  0);  (7) 

whenever  P/a(&}( Y)  |  0  =  0)  <  a,  9  denotes  the  actual  frequency  of  the  signal.*  The  rule  6?(Y)  is 
assumed  to  be  of  the  form 


1 


MY)  = 


<  0  or  1  if  max,  |  CjY 


0 


> 

=  r, 
< 


(8) 


•We  assume  Pd,  Pfa  can  be  written  in  terms  of  7  to  not  depend  on  s  or  a  separately. 
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where  i  =  1 , . . .  ,r,  r  being  the  number  of  filters.  Each  filter  coefficient  vector  C,  will  be  associated 
with  a  frequency  range  of  interest  A,  C  [0,  2i r]  such  that  (A(,  . . Arj  forms  a  partition  of  the  inter¬ 
val  [0.  lit],  Suppose  a  filter  with  coefficients  C  is  to  be  designed  for  and  interval  A^  €  (Aj, Ar|. 
It  is  desired  that  Pd(5?( Y)  |  0  =  8)  be  close  in  some  sense  to  the  optimum  Pd(8}( Y)  |  0  =  8)  for  each 
B*\a. 

Definition  1  —  Define  yc  as  the  input  SIR  required  such  that  5?‘(Y)  will  have  the  same 
Pq  as  5J(Y)  at  0  =  8,  i.e., 

PMC( Y)  |  0  =  8)  =  Pdm Y)  |  0  =  0), 

when  Pfa  =  a  for  both  tests.  The  log  energy  e(6)  is  given  as 

e(6)  =  log  yc  -  log  y. 

Furthermore  the  log  energy  risk  is  defined  as 


r(8)  =  max  g(8)e  (8), 

0«a„ 

where  g{8)  is  an  arbitrary  weighting  function. 

The  log  energy  function  e(8)  is  the  log  energy  required  for  the  filter  with  coefficients  C  to  per¬ 
form  the  same  as  the  optimal  filter  at  the  given  frequency  8.  The  log  energy  risk  is  defined  as  the 
minimum  log  energy  needed  for  the  suboptimal  detector  to  perform  at  least  as  well  as  the  6$  (Y)  for 
all  6  6  A,,.  This  risk  can  be  weighted  depending  on  the  cost  structure  desired  on  the  set  Aa.  The 
weighting  will  take  the  form  of  a  function  g(0).  The  risk  can  be  minimized  in  a  Bayesian  sense, 

5  €  arg  min  max  g(8)e(8) 
s  ei\A 

or  this  can  be  written  in  terms  of  the  filter  C, 

C  €  arg  min  maxg(0)e(0).  (9) 

c  Si\A 

The  log  energy  e(8)  accurately  reflects  the  degradation  of  the  desired  filters  relative  to  what  can 
be  achieved  in  the  optimal  case  in  a  manner  relevant  to  radar  engineering.  The  number  of  filters  that 
can  be  designed  by  using  this  criterion  is  not  constrained,  nor  is  the  frequency  separation  of  the 
filters,  as  the  intervals  A  can  arbitrarily  be  chosen.  Thus,  the  criterion  is  broad  enough  to  allow 
important  design  parameters  to  be  varied.  Finally,  the  criterion  is  capable  of  being  weighted  by  the 
arbitrary  function  g(8),  which  allows  it  to  incorporate  desired  cost  functions  such  as  velocity  profiles. 

It  is  interesting  to  compare  this  criterion  with  the  maximum  formulation  used  in  many  robustness 
problems  (12,13]  in  which  one  seeks  a  filter  to  maximize  the  minimum  gain  over  a  class  of  possible 
signals  and/or  interference  models. 


max  min  SIR .  (10) 

c  ei\A 

This  theory  has  been  developed  so  that  explicit  solutions  are  available  under  certain  hypothesis.  In 
our  case  the  signal  may  vary  by  taking  on  different  6.  This  criterion  is  compared  with  the  MMLE  ir. 
the  Section  5  for  the  design  of  single  filters. 
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4.  MMLE  APPLIED  TO  THE  RADAR  PROBLEM 


To  apply  the  MMLE  to  the  radar  problem  Eq.  (1),  e(6)  must  be  determined  in  Eq.  (9).  This 
requires  the  rules  5j)(Y)  and  the  associated  Pd  and  V  The  most  powerfull  tests  6}  by  the  Neyman- 
Pearson  criterion  are  found  to  be  a  form  of  the  Wiener  filter. 


1 


5(Y)  = 


0  or 
0 


> 

if|YrR_1S*|  =  7(6). 

< 


01) 


The  Pd  for  this  test  is  Ref.  14 


where  /  =  |  CTS  \  ,  r2  =  CfRC ,  Q  is  the  Marcum  Q  function.  We  wish  to  find  the  input  SIR  yc 
such  that  the  suboptimal  rule  6^(Y)  performs  the  same  as  the  test  5^(Y)  for  some  6  €  Aa.  This 
requires, 


Q 


=  Q 


6  €  , 


where  lc,  tc  correspond  to  the  test  yy/  (Y)  and  1(6),  t(9)  to  the  respective  b}(\).  This  amounts  to 


K6) 

e(6) 


lc 

—  6  £  Aa, 

Tc 


which  is  equivalent  to 


where 


Gel  ~  Gdc . 


i  1 2 

G(6)  =  -  \  1  (12) 

7  C^C 

is  the  gain  for  a  filter  with  coefficients  C,  and  +  denotes  conjugate  transpose.  Ge,  Gc  denotes  the 
gains  for  the  most  powerful  test  and  the  suboptimal  test  respectively.  Now  e(6)  can  be  determined 
from 

e(6)  =  log  yc  -  log  7  =  log  Ge  -  log  Gc, 
and  substituting  for  the  coefficients  of  the  optimal  test  6^(Y)  yields 


=  log  StR-'S*  Cr/?C2 

e  |  crs  1 2 


(13) 


6 


NRL  REPORT  9280 


To  find  the  log  energy  risk,  e(6)  is  maximized  with  respect  to  9.  Assuming  uniform  costs,  i.e., 
g(6)  =  1  V0  €  A^, 


r(S)  =  max  e(6) 
ot\A 


=  C  RC  max 

6€Aa 


StR~1S* 

I  CrS  1 2  ' 


(14) 


To  find  the  rule  or  suboptimal  coefficient  vector  C,  the  log  energy  risk  is  minimized. 

Sr/?  _1S* 


min 

c 


Cf/?C  max  T  , 
BtAA  C7S  I  2 


(15) 


A  filter  with  coefficients  C  is  obtained  when  the  minimum  is  reached  (it  may  or  may  not  be 
unique).  We  call  filters  derived  in  such  a  manner  MMLE  filters.  Note  that  the  maximum  over  6  is 
essentially  the  maximum  of  the  ratio  of  two  polynomials  with  coefficients  from  C  or  R.  No  general 
solution  appears  other  than  a  search  algorithm.  Numerical  algorithms  have  been  developed  and  are 
discussed  in  Appendix  C. 

5.  DESIGN  FOR  LOW-DIMENSION  N 

The  results  shown  in  the  preceding  section  will  now  be  used  to  design  a  detector  incorporating  a 
single  filter,  when  the  dimension  N  is  low.  The  degree  to  which  N  is  considered  low  is  relative  to  the 
particular  interference  covariance  matrix  under  consideration,  hence  may  change  for  different  applica¬ 
tions. 


We  will  consider  three  scenarios  and  the  corresponding  covariance  matrices. 

(a)  The  received  interference  is  dominated  by  clutter  introduced  by  rain,  cloud,  foliage,  and  the 
like.  The  clutter  is  highly  correlated  and  will  be  modeled  as  a  Gaussian  process  with  a  Gaussian- 
shaped  clutter  power  spectrum.  The  covariance  matrix  can  be  written  as  in  Ref.  15, 

1  +  a2x  p(T)  ■■■  p(NT) 

p(T)  l  +  o~x  ’ ' '  : 

p(NT)  p((N  -  1) T)  ■■■  1  +  a2x 

where 

p(kD  =  e'2"k0cT)' .  (17) 

The  thermal  noise  energy  per  dimension  is  a\ ,  the  clutter  energy  per  dimension  is  normalized  to 
1,  and  the  pulse  repetition  is  T.  The  total  interference  energy  per  dimension  a2  is  then  1  +  a2x.  The 
standard  deviation  of  the  Gaussian  clutter  power  spectrum  is  ac ,  and  the  normalized  spectral  width  of 
the  spectrum  is  oqT. 
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(b)  Extraneous  interference  in  a  radar  return  may  cause  one  or  more  components  to  be  cor¬ 
rupted.  For  simplicity,  one  component  will  be  assumed  corrupted.  It  is  assumed  that  the  interference 
dominates  so  that  no  useful  information  can  be  extracted  from  the  corrupted  pulses  in  the  detection 
process.  The  optimal  solution  for  a  fixed  frequency  and  /'th  pulse  corrupted  is  to  remove  the  /th  row 
and  /'th  column  of  the  covariance  matrix  above  and  accordingly  process  the  remaining  pulses  [16]. 
The  operator  Pj(X )  denotes  the  vector  X  with  the  /'th  component  removed. 

c.  A  covariance  matrix  is  devised  in  Appendix  A  to  illustrate  the  Emerson  procedure.  This 
matrix  can  be  written  as, 


S,S,+ 

*  =  Ex<i r  s'  = 


.  2ir(i  -  1 ) 
N 

,  4ir(i-i) 
N 


.  2t|N -!)(/- 1) 
N 


(18) 


where  /  =  1  . . .  ,  N.  The  set  [X] ,  . . .  ,  X^l  contains  the  eigenvalues  of  the  matrix  R,  and  the  eigen¬ 
values  are  assumed  distinct.  Note  that  the  matrix  is  Toeplitz. 


For  each  of  the  three  cases,  the  Emerson  and  the  MMLE  methods  are  compared  with  the  GLRT 
or  infinite  filter  bank,  Eq.  (6).  The  gain  function  for  the  GLRT  is  optimal  for  a  given  6.  Henceforth 
if  is  referred  to  as  the  optimal  gain.  The  GLRT  as  a  detector  however  is  not  optimal  when  one  does 
not  have  a  priori  knowledge  of  6.  We  assume  that  the  speed  of  the  target  is  unknown  and  is  uniformly 
distributed  in  the  interval  [0,  2i r].  Furthermore  we  assume  a  uniform  cost  function  g(8)  in  Eq.  (9). 
For  simplicity,  oqT  is  fixed  at  0.05  and  C/N0  at  50  dB.  The  design  goal  is  to  be  as  close  to  the 
optimal  gain  in  the  passband  as  possible.  Considering  acT  =  0.05,  a  somewhat  arbitrary  passband 
will  be  assumed  between  the  normalized  frequency  (f  =  d/2ir)  intervals  /  €  [1,  .9]  and  /  €  [.2,  .8] 
depending  on  the  remaining  parameters.  If  a  single  filter  is  not  sufficient  to  cover  the  passband  satis¬ 
factorily,  it  is  necessary  to  either  change  the  passband  or  add  additional  filters  as  discussed  in  the  next 
section. 


In  Fig.  1,  gain  is  plotted  vs  /,  for  the  covariance  matrix,  Eq.  (16)  for  N  =  4.  The  optimal  gain 
obtained  by  using  the  filter  coefficients  R~l S*  is  plotted  along  with  the  single  filter  gains  of  Eq.  (12) 
based  on  Emerson's  weights  and  the  MMLE.  It  is  clear,  in  the  passband,  that  both  single  filters  per¬ 
form  nearly  the  same  as  the  GLRT,  Eq.  (6).  Nearly  the  same  performance  is  observed  for  N  =  6 
shown  in  Fig.  2.  except  the  MMLE  begins  to  separate  from  Emerson's  filter  at  the  beginning  and 
middle  of  the  passband.  This  separation  can  be  seen  clearly  for  N  =  8  in  Fig.  3  where  the  MMLE 
achieves  at  most  20  dB  more  improvement  over  Emerson's  method  for  frequencies  in  the  intervals 
/  6  [.15,  .35],  /  €  [.65,  .85],  and  sacrifices  at  most  5  dB  in  the  interval  /  €  [.35,  .65].  The  maxi¬ 
min  formulation  of  Eq.  (10)  is  shown  in  Fig.  4  relative  to  the  MMLE  and  infinite  bank  Eq.  (6)  for 
N  =  8.  Both  the  MMLE  and  maximin  were  optimized  over  /  6  [  2,  .8],  The  relatively  poor  perfor¬ 
mance  of  the  maximin  is  heuristically  due  to  maximizing  the  minimum  gain  that  occurs  at  /  =  .2  and 
/  =  .8  in  the  interval  /  €  [.2,  .8],  This  is  done  at  the  expense  of  the  remaining  passband  that  is 
allowed  to  sharply  degrade.  A  more  relevant  application  of  the  maximin  procedure  is  discussed  in 
Section  8. 
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-  OPTIMAL 

- MULE 

-  EMERSON 


o’.l  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9 


NORMRLIZED  FREQUENCY 

Fig.  1  —  Gain-Gaussian  clutter  power  spectrum  N  =  4. 
acT  =  0.05.  C/N0  =  50  dB 


OPTIMAL 

MMLE 

EMERSON 


Fig.  2  —  Gain-Gaussian  clutter  power  spectrum  N  =  6. 
acT  =  0.05,  C/N0  =  50  dB 
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Fig.  3  —  Gain-Gaussian  clutter  power  spectrum  N  =  8. 
ocT  =  0.05.  C/N0  =  50  dB 
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In  Fig.  5,  curves  of  the  optimal  gain,  MMLE  gain,  and  gain  obtained  by  using  Emerson’s  Filter 
are  shown  for  N  =  5.  The  third  row  and  column  of  Eq.  (16)  are  removed  to  account  for  extraneous 
interference.  The  gain  function  for  Emerson's  filter  yields  a  blind  speed  at  /  =  .5,  which  could  be 
critical  for  many  systems  (this  problem  has  been  observed  for  odd  M).  The  MMLE  maintains  a 
nearly  identical  loss  across  the  passband.  The  single  MMLE  filter  again  yields  uniform  coverage  in 
Fig.  6  for  N  =  9  while  Emerson's  filter  produces  poor  results  at  /  =  0.5. 


Fig.  5  —  Gain-extraneous  interference  Py.  N  =  5,  acT  =  0.05.  C/A/n  =  50  dB 
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We  now  consider  the  covariance  matrix  of  the  form  of  Eq.  (18).  Again  Emerson's  method  and 
the  MMLE  are  compared  with  the  infinite  bank.  For  N  =  3  shown  in  Fig.  7,  the  eigenvalues  are 
chosen  so  that  most  of  the  interference  is  near  zero  Doppler;  X,  =  1.,  X2  =  0.01,  and  X3  =  0.0099. 
As  explained  in  Appendix  A,  choosing  any  eigenvector  of  the  covariance  matrix  results  in  two  nulls 
in  the  Doppler  space,  in  our  case  at  /  =  0,  /  =  2/3.  The  gain  for  the  MMLE  design  is  uniformly 
close  to  the  optimal  gain.  The  same  effect  is  seen  in  Fig.  8  in  which  N  =  5.  Here  X]  =  1., 
X2  =  0.0000999,  X3  =  0.0001,  X4  =  0.0001,  X5  =  0.0001.  These  cases  are  based  on  a  covariance 
matrix  that  is  designed  to  help  in  understanding  the  problem  and  may  not  be  indicative  of  covariance 
matrices  encountered  in  physical  processes. 


Fig.  7  —  Gain-interference  matrix  (17).  N  =  3 


Fig.  8  —  Gain-interference  matrix  (17),  (V  =  5 
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6.  DESIGN  FOR  HIGH-DIMENSION  N 

When  the  single  filter  designs  of  the  previous  section  fail  to  achieve  the  desired  goal,  a  useful 
recourse  is  to  consider  the  implementation  of  multiple  filters.  To  compare  the  designs,  we  define  the 
coverage  C(/)  as. 


C(f)  =  max  G,(2irf),  i  =  1 . P  (19) 

/ 

where  G,(27r/)  is  the  associated  gain  function,  Eq.  (12)  of  the  ith  filter,  and  P  is  the  number  of 
filters.  The  uniform  weighting  function,  g(6)  =  1  in  the  passband;  g(6)  =  0  otherwise,  is  assumed. 
It  is  desirable  to  have  the  coverage  as  close  in  the  passband  to  what  is  obtainable  with  the  infinite 
bank  of  filters,  Eq.  (6)  where  each  filter  is  optimal  for  a  corresponding  frequency.  The  covariance 
matrices  that  are  considered  are  cases  1  and  2  of  Section  5,  that  of  the  interference  having  a  Gaussian 
power  spectrum  and  extraneous  interference  corrupting  one  component  of  the  received  vector.  Again 
we  assume  acT  =  0.05  and  C /A/0  =  50  dB. 

Comparisons  are  made  between  the  MMLE  the  N  filters  by  Andrews,  and  the  approximated 
CGLRT  with  N  equally  spaced  filters.  Andrews’  design  for  multiple  filters  is  a  straightforward  exten¬ 
sion  of  the  Emerson’s  procedure.  Here  the  Doppler  space  is  partitioned  into  N  equally  spaced  inter¬ 
vals,  and  N  filters  are  found  by  maximizing  the  average  gain  for  each  respective  interval.  Specifically 
the  ith  interval  is  given  by 


"(2  i  -  1) 
2  N 


(2/  +  1) 

2  N 

J 


(20) 


where  /  =  0 . N  -  1  and  the  units  are  normalized  frequency.  This  is  in  contrast  to  Emerson’s 

procedure  that  yields  a  single  filter  by  maximizing  the  average  gain  over  the  entire  [0,  1]  normalized 
frequency  interval.  The  CGLRT  consists  of  N  filters  with  coefficients  of  the  form  By  parti¬ 

tioning  the  Doppler  space  into  N  equally  spaced  intervals  in  Eq.  (20)  and  by  using  the  average  fre¬ 
quency  in  each  interval  to  evaluate  the  signal,  N  such  filters  are  derived.  The  MMLE  is  well  suited 
for  the  problem  of  multiple  filter  design,  since  filters  are  optimized  over  disjoint  frequency  intervals 
A,  as  described  in  Section  3.  To  meet  a  design  goal  with  the  least  number  of  filters,  one  must  deter¬ 
mine  the  frequency  intervals  A,  .  Given  the  same  general  design  goals  and  passbands  as  in  Section  3, 
this  can  be  done  as  follows.  Assuming  a  single  filter  did  not  meet  the  design  requirements,  the  first 
interval  A|  is  set  equal  to  the  interval  /  €  [0,  .5].  The  MMLE  is  applied  and  the  gain  is  determined. 
If  the  gain  is  not  sufficiently  close  to  the  optimal  gain  in  the  interval  /  6  [0,  .5],  the  interval  is 
decreased  to  /  €  [0,  jc]  where  x  is  less  then  0.5,  and  the  MMLE  is  applied  again.  Changes  (both 
increasing  and  decreasing  as  necessary)  are  continued  in  this  manner  until  the  gain  meets  the  design 
requirements  in  the  interval.  Suppose  the  single  filter  now  meets  the  requirements  in  the  interval 
/  €  [0,  jc],  and  a  second  filter  is  designed  in  the  interval  /  €  [jc,  1]  and  the  coverage,  Eq.  (19),  of  the 
first  and  second  filters  is  compared  with  the  design  requirements.  This  process  of  changing  the  inter¬ 
val  is  continued  as  above  (/  €  [jc,  y  ]  y  variable)  until  the  coverage  is  satisfactory  through  the  second 
interval.  Enough  filters  are  added  by  following  this  procedure  until  the  passband  is  covered  to  meet 
the  requirements 


Figure  9  shows  the  approximated  CGLRT  vs  the  GLRT  or  infinite  filter  bank,  Eq.  (6),  for 
A  =  10.  The  coverage  provided  in  the  passband  is  practically  identical  to  the  GLRT  coverage.  This 
coverage  is  repeated  in  Fig.  10  in  which  filters  are  designed  by  using  the  Andrews’  method.  Finally 
the  MMLE  is  shown  in  Fig.  11  in  which  only  four  filters  are  used.  The  intervals  determined  by  the 
method  above  are  /  €  (.17,  .4]  ,  /€  [  4,  .5);  and  conjugate  filters  are  used  for  the  intervals 
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o  c.2  o'.:  c.5  0.5  0.7  0.8  o'. 3  i 


NORMAL :iEw  EREQUENCY 

Fig.  9  —  Coverage-Gaussian  clutter  power  spectrum  multiple  filters  CGLRT. 
N  =  10,  ocT  =  0.05,  C/N„  =  50  dB 


a 


Fig.  10  —  Coverage-Gaussian  clutter  power  spectrum  multiple  filters  Andrews, 
N  =  10.  ocT  =  0.05,  C/A',,  =  50  dB 
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o 


normalized  frequency 

Fig.  1 1  —  Coverage-Gaussian  clutter  power  spectrum  multiple  filters  MMLE, 
N  =  10,  acT  =  0.05,  C/N0  =  50  dB 


/  €  [.5,. 6]  ,  /  €  [.6, .83].  Here  the  length  of  the  interval  of  the  first  filter  is  0.23  while  the  second  is 
less  then  half  of  this  at  0.1.  Nonuniform  lengths  are  important  if  a  minimum  number  of  filters  is 
desired. 

Extraneous  interference  is  considered.  Figures  12  and  13  compare  the  CGLRT  and  Andrews’ 
filters  respectively  to  the  infinite  bank.  The  dimension  N  =  9  and  the  middle  component  is  corrupted, 
resulting  in  eight  pulses.  The  coverage  provided  in  the  passband  is  nearly  the  same  as  the  infinite 
filter  bank  coverage.  Figure  14  shows  the  MMLE  design  in  which  only  three  filters  were  needed  to 
approximate  the  optimal  gain.  The  intervals  used  in  the  design  were  /  6  [17,  .38],  its  conjugate  filter 
interval  /  €  [.62,  .83],  and  /  €  [.38,  .62].  The  lengths  of  the  intervals  are  nearly  equal. 


Fig.  12  —  Coverage-extraneous  interference  P5,  multiple  filters  CGLRT, 
N  =  9,  acT  =  0.05,  C/Na  =  50  dB 
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Fig.  13  —  Coverage-extraneous  interference  P5.  multiple  filters  Andrews. 
N  =  9,  ocT  =  0.05.  C/N0  =  50  dB 


Fig.  14  —  Coverage-extraneous  interference  P multiple  filters  MMLE. 
N  =  9,  trcT  =  0.05.  C/Nn  =  50  dB 
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In  summary  we  have  shown  that  only  four  filters  are  needed  to  effectively  cover  the  p^ssband 
for  N  =  10.  For  extraneous  interference  only  three  filters  are  necessary. 

7.  ROBUSTNESS 

The  filters  in  the  previous  section  were  designed  under  the  assumption  that  the  statistics  of  the 
interference  are  accurately  reflected  by  the  interference  covariance  matrix  used.  This  may  be  so  in 
certain  space  based  applications.  However,  for  other  applications  this  may  not  be  so.  In  fact,  the 
interference  matrix  may  change  from  range  cell  to  range  cell  and  from  scan  to  scan  of  a  pulsed  radar. 
The  received  signal  may  either  possess  Gaussian  or  nongaussian  statistics.  For  our  purposes  a  Gaus¬ 
sian  process  is  assumed  and  furthermore  we  assume  that  the  received  signal,  Eq.  (1)  (which  is  at  a 
constant  range),  is  wide  sense  stationary.  The  actual  covariance  matrix  associated  with  the  Gaussian 
interference  is  assumed  to  lie  within  a  class  of  covariance  matrices  that  have  a  Gaussian  clutter  power 
spectrum.  The  class  3  of  covariance  matrices  is  defined  as 


I?  as  in  Eq.  (16)  |  ocT  €  Aa 


(21) 


where  Aa  is  the  set  of  relevant  normalized  spectral  widths,  and  Ah  is  the  set  of  relevant  clutter  to 
thermal  noise  ratios,  C/N0  =  \/o\.  Relationships  on  filter  performance  when  C/N0  is  varied  are 
presented.  The  robustness  of  the  MMLE,  Andrews'  method,  and  the  CGLRT  filters  will  be  examined 
when  they  are  designed  at  a  fixed  operating  point  of  ocT  and  C/N o- 

We  begin  with  an  equation  derived  from  Ref.  7  that  relates  gain  averaged  over  the  entire  nor¬ 
malized  frequency  interval  [0,  1],  termed  the  improvement  factor  IF,  when  Emerson’s  method  is 
used. 


_1_ 

len 


c 

N0 


1  + 


C 

Nn 


1 


1  + 


C 

No 


(22) 


/(.  represents  the  IF  when  the  eigenvector  corresponding  to  the  minimum  eigenvalue  of  the 
clutter  covariance  matrix  is  used  as  the  filter  weight  ( C/N0  =  oo).  By  use  of  the  same  filter,  lcn  is 
the  new  IF  when  white  noise  is  added  to  the  clutter.  By  taking  logarithms  of  Eq.  (22)  and  approxi¬ 
mating  for  C/Nq  »  1, 


10  log  Icn 


10  log  — - 10  log 

No 


f — 
No 

lc 


+  1 


(23) 


If  C/N o  «  /f,  then  10  log  /r„  =  10  log  C/N0.  If  C/N0  »  7f,  then  10  log  lcn  ~  10  log 
Ic.  Hence  C/N o  is  a  limiting  factor  in  the  improvement  factor  when  white  noise  is  added.  Although 
this  relation  yields  a  useful  result,  it  applies  to  Emerson  filters  and  does  not  hold  for  general  linear 
filters.  Furthermore,  Eq.  (22)  does  not  tell  us  what  happens  as  a  function  of  frequency,  only  the 
average  gain  or  improvement  factor  is  presented.  A  more  general  result  is  stated  here  as  a  robustness 
theorem. 
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Theorem  1  —  Given  a  vector  of  filter  coefficients  let  Gc(6 )  represent  the  gain,  Eq.  (12)  for  which  R  is 
a  N  X  N  positive  definite  clutter  covariance  marix  and  diagonal  elements  a].  Assuming  the  same 
filter  coefficients,  let  Gcn(6)  represent  the  gain  for  which  the  covariance  matrix  is  R'  =  R  + 

yg|  -T-  1.  IfA(d)  =  10  log  Gc(6)  -  10  log  Gcn(6), 

C 


10  log 


C 

No 


1  + 


No 


1  + 


where  Gmax  =  maxsGr(0). 


N 


Nr 


<  A (6)  <  10  log 


C 

N0 


1  + 


Nn 


1  + 


C 

No 


(24) 


For  a  proof  see  Appendix  B.  If  C/N0  >  Gmax.  the  upper  bound  indicates  A (0)  is  small.  So 
long  as  C/Nq  >  G^,  any  linear  filter  is  robust  in  the  sense  that  the  gain  function  will  not  change 
as  a  function  of  C/N0.  If  additionally  Gc(d)  is  close  to  the  optimal  gain  in  which  C/N0  =  oo,  then 
Gcn(0)  will  also  be  close  to  the  optimal  gain  for  every  C/N0  >  Gmax.  This  follows,  since  the 
optimal  gain  for  any  R'  with  finite  C/N0  is  less  then  the  optimal  gain  with  covariance  matrix  R.  In 
the  case  where  C/N0  <  Gmax,  Eq.  (24)  yields  upper  and  lower  bounds  on  the  degradation  and  the 
gain  function  will  suffer.  These  bounds  are  plotted  in  Fig.  15  for  N  =  5,  acT  =  0.05,  and 
C/N o  =  10,  30,  oo  dB.  The  topmost  figure  is  Gc  for  a  single  filter  designed  by  the  MMLE  for 
C/Nq  =  oo.  The  curves  portrayed  with  C/N0  =  30  dB  represent  the  same  filter  but  with  dif¬ 

ferent  covariance  matrices  incorporating  the  respective  white  noise.  The  dotted  line  curves  above  and 
below  the  C/N0  =  10,  30  dB  curves  are  the  upper  and  lower  bounds  in  Eq.  (24).  The  bounds  hold 
whether  or  not  the  interference  is  Gaussian. 


Fig.  15  —  Robustness  bounds,  S  =  5,  ocT  +  0.05,  C/N0  -  10.  30.  oo  dB 
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Now  we  address  the  question  of  robustness  for  the  filters  considered  in  the  previous  sections. 
For  comparison,  the  filters  are  all  designed  assuming  acT  =  0.05,  C/N{)  =  90  dB,  with  passband 
/  €  (.2.  .8]  and  uniform  weighting  function  in  the  passband.  The  dimensionality  will  be  fixed  at 
N  =  10.  The  degradation  of  each  filter  relative  to  the  infinite  bank  of  filters  for  different  o^T,  C/N 0, 
and  /.  will  be  shown. 

Figure  16  is  a  robustness  curve  of  a  bank  of  ten  CGLRT  filters  designed  as  discussed  in  Section 
6  at  /  =  0.5.  Each  point  of  the  solid  lines  represents  an  optimal  filter  with  weighting  R~lS* 
matched  at  the  particular  acT,  C/N0  and  /.  The  dotted  lines  are  the  coverage  that  is  obtained  with 
the  fixed  bank  of  CGLRT  filters  designed  with  the  parameters  above.  This  filter  bank  is  particularly 
robust  for  C/N0  <  80  dB  and  the  large  range  of  acT  shown.  On  the  other  hand,  when  /  =  .2,  Fig. 
17,  the  results  are  rather  poor  relative  to  what  can  theoretically  be  achieved.  In  fact,  there  are  regions 
where  the  coverage  could  increase  by  30  dB  or  more  by  better  processing.  Nearly  the  same  coverage 
is  seen  in  Figs.  18  and  19  for  /  =  0.5,  /  =  0.2  respectively  when  the  design  criterion  is  Andrews’ 
extension  of  Emerson's  filters.  Only  four  filters  are  used  in  the  MMLE  design  as  shown  in  Fig.  20 
for  /  =  0.5.  The  filters  were  designed  by  partitioning  the  passband  as  /  €  |.2,  .25],  /  6  ] .25,  .5], 
/  €  [,5,  .75],  /  €  [.75,  .8],  The  filters  are  particularly  robust  at  /  =  0.5  for  the  parameters  shown. 
In  Fig.  21  the  coverage  is  shown  for  /  =  0.2.  As  with  the  previous  cases,  the  coverage  is  poor  com¬ 
pared  with  the  optimal  gain. 

The  class  of  filters  known  as  Chebychev  is  designed  somewhat  independent  of  ocT,C/N0 ,  and 
/.  The  Chebychev  filters  can  be  characterized  by  their  sidelobe  levels  and  dimensionality  [17],  In 
Fig.  22  the  coverage  is  shown  for  ten  Chebychev  filters,  equally  spaced,  Eq.  (20).  at  f  =  0.5  with  90 
dB  sidelobe  levels.  The  coverage  is  mostly  identical  to  the  coverage  obtained  with  the  previous  filters. 
When  /  =  0.2  as  in  Fig.  23,  the  coverage  is  much  better  than  that  obtained  in  the  previous  cases. 
Unfortunately  it  still  lacks  the  optimal  coverage  by  at  most  7  dB  at  C/N0  =  20  dB,  8  dB  at 
C/N0  =  50  dB,  and  23  dB  at  C/N0  =  80  dB. 

CD 


Fig.  16  -  Coverage  CGLRT  /  =  0.5.  N  =  10  a(  T  =  0.05. 
C/N„  =  20.  50.  80.  <x  dB 
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Fig.  21  —  Coverage  MMLE  /  =  0.2.  N  =  10  acT  =  0.05. 
C/N0  =  20.  50.  80.  oo  dB 


a  i 


Fig.  22  —  Coverage  Cheby  /  =  0.5,  N  -  10 
CfNo  =  20.  50.  80,  oo  dB 
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Fig.  23  —  Coverage  Cheby  /  =  0.2.  N  =  10 
C/N0  =  20.  50,  80,  oo  dB 


The  problem  of  obtaining  filters  with  frequency  coverage  nearly  optimal  for  a  particular  operat¬ 
ing  point  of  C/N0  and  acT  has  been  solved.  When  C/Nq  and  are  allowed  to  vary,  such  filters 
have  been  shown  to  degrade  relative  to  the  optimal  coverage.  This  is  further  examined  in  the  next 
section. 

8.  CONCLUSIONS 

We  have  presented  a  design  method  that  will  produce  filters  capable  of  detecting  targets  of 
unknown  velocity.  A  design  criterion  is  specified  by  the  radar  designer  in  terms  of  gain  relative  to 
the  optimal  gain  as  a  function  of  frequency.  The  method  requires  partitioning  of  the  passband  and 
comparison  of  the  coverage  of  the  corresponding  MMLE  filters  with  the  design  criteria.  Enough 
filters  are  added  as  described  in  Section  6  until  the  coverage  meets  the  design  criteria.  The  iterative 
partitioning  is  done  in  an  attempt  to  minimize  the  number  of  filters  needed.  The  degree  of  the  close¬ 
ness  was  explained  in  Section  3.  The  iterative  partitioning  will  very  likely  produce  nonuniform  length 
intervals. 

These  filters  were  designed  under  the  premise  of  known  values  of  acT  and  C/N0.  The  case, 
when  these  values  are  not  known  precisely,  is  discussed  in  Section  7.  General  bounds  were  presented 
for  a  given  linear  filter  when  the  clutter  to  thermal  noise  ratio  is  allowed  to  vary.  The  performance  of 
filters  such  as  the  MMLE,  which  require  precise  knowledge  of  the  covariance  matrix,  was  deter¬ 
mined.  This  was  addressed  by  introducting  a  class  of  covariance  matrices  in  which  ocT  and 
C/Nq  are  allowed  to  vary.  The  actual  design  of  filters  that  are  robust  within  the  covariance  matrix 
class  ^  is  a  remaining  problem.  An  important  step  in  designing  such  filters  is  to  define  a  criterion 
appropriate  to  this  problem.  Obviously  the  usefulness  of  the  maximin  formulations, 

max  min  SIR 
c  rt 
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where  the  maximization  is  over  the  filter  coefficients  C.  and  MMLE  needs  to  be  determined.  Addi¬ 
tionally  alternative  criteria  [15]  should  be  investigated. 
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Appendix  A 

IMPROVEMENT  FACTOR 


The  criteria  of  maximizing  the  IF  (improvement  factor),  referred  to  as  Emerson’s  criterion,  is 
now  examined.  The  primary  consideration  is  to  determine  when  IF  is  a  good  indicator  of  filter  per¬ 
formance.  To  better  understand  the  criteria,  the  covariance  matrix. 


f  1 

.  2t(i  - 1 ) 


s,s; 


'  N 
.  4t (/  -  1 ) 
1  N 


.  2ir{N  —  !)(/  —  1 ) 
7 - N - 


where  /  =  1 . N  will  be  applied.  The  set  {X, . \N\  contains  the  eigenvalues  of  the  matrix. 

The  eigenvalues  arc  assumed  distinct  and  the  matrix  is  Toeplitz.  We  set  the  eigenvalues  to 
=  1.,  X2  =  0.00001,  X-,  =  0.0000099,  X4  =  0.00001,  X5  =  0.00001.  Note  that  when  Gaussian 
interference  associated  with  this  covariance  matrix  is  passed  through  a  linear  filter,  the  IF  is  not 
1/Xmin,  where  Xmin  is  the  minimum  eigenvalue.  This  is  because  the  diagonal  elements  of  R  are  not 

necessarily  equal  to  one.  Assuming  the  same  hypotheses  as  in  Eq.  (1)  IF  = - ,  where  7  is  the 

7  7min 

input  SIR.  The  solid  curve  in  Fig.  24  is  the  gain  Eq.  (12)  of  the  infinite  filter  bank,  Eq.  (6).  Also 
shown  are  curves  obtained  through  the  MMLE  and  Emerson's  criteria.  Emerson's  criterion  yields  a 
filter  with  rather  poor  gain  response.  This  is  because  the  eigenvector  corresponding  to  the  minimum 

eigenvalue  must  be  S,  for  some  /.  Every  eigenvector  is  orthogonal  to  every  other  eigenvector. 

Hence,  the  numerator  of  Eq.  (12)  is  zero  whenever  the  conjugate  of  S  is  orthogonal  to  the  eigenvec¬ 
tor.  This  happens  N  -  1  times  in  our  example.  Essentially  CrS  where  6  €  [0,2ir]  represents  the  con¬ 
jugate  of  the  discrete  time  Fourier  transform  [Al]  of  the  coefficient  vector  C.  Consider  the  z 
transform  LjCjZ~j ,  where  c  is  a  complex  variable.  This  equation  has  N  -  1  roots 

Zj,  i  =  1,  2 . N  —  1.  If  a  root  is  on  the  unit  circle  in  the  z  plane  z,  =  e7(3,  the  gain  will  be 

forced  to  zero  at  the  frequency  6  =  — /3.  Robinson  [A2]  showed  that  all  the  roots  of  the  eigenvector 
corresponding  to  the  minimum  eigenvalue  are  on  the  unit  circle  when  the  covariance  matrix  is  Toe¬ 
plitz.  This  may  be  useful  if  the  roots  are  close  to  6  =  0  since  it  will  notch  the  clutter.  However,  in 

cases  such  as  in  Fig.  Al,  much  of  the  remaining  Doppler  space  is  also  notched. 

Emerson’s  criterion  consists  of  finding  the  normalized  filter  coefficients  that  minimize  the 
interference  output  power.  Equivalently  the  average  gain  or  IF  can  be  maximized.  In  the  derivation 
of  Emerson’s  filter,  we  assume  that  the  Doppler  speed  is  uniformly  distributed.  Also,  no  weighting  or 
cost  function  of  frequency  exists,  i.e.,  a  uniform  cost  is  assumed.  The  coarse  design  goal  is  to  pro¬ 
duce  a  filter  that  delivers  favorable  gain  coverage  under  these  circumstances.  One  measure  for  judg¬ 
ing  the  effectiveness  of  such  a  criterion  is  to  determine  how  the  IF  relates  to  the  design  goal.  Since 
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the  IF  is  to  be  maximized,  does  a  higher  IF,  in  fact,  indicate  a  better  design?  Figure  A1  shows  that 
the  MMLE  yields  a  filter  design  uniformly  close  to  the  optimal  gain.  It  outperforms  Emerson's  filter 
in  80%  of  the  Doppler  space  and  does  not  suffer  the  numerous  holes  in  the  coverage.  Yet,  for  the 

MMLE  design,  IF  =  42.997  dB,  and  for  Emerson’s  design  IF  =  43.000  dB.  In  this  example,  the  IF 

does  not  effectively  relate  to  the  design  goal.  To  further  illustrate  this,  examine  Fig.  A2  forN  =  5, 
acT  =  0.05,  c/N0  =  50  dB.  The  solid  curve  was  obtained  from  a  single  filter  design  by  MMLE 
with  passband  /  6  [.2,. 8]  and  uniform  weighting.  It  was  determined  the  IF  =  43.76  dB.  Also  shown 
in  the  figure  by  the  dotted  lines  is  a  rectangular  window.  This  window  has  no  gain  in 
/  6  [0,  .4],  /€  [.6,1]  and  performs  slightly  better  then  the  MMLE  in  /  6  [.4,  .6],  Obviously  this 

would  make  a  rather  poor  filter  compared  with  the  MMLE.  The  IF  is  43.82  dB,  which  is  slightly 

higher  then  the  MMLE  IF.  Although  the  actual  eigenvector  solution  performs  much  better,  and  an 
actual  filter  that  yields  the  dotted  gain  curve  may  not  even  exist,  the  example  does  help  to  show  that 
IF  is  not  necessarily  a  good  indicator  of  the  design  goal. 

In  summary,  Emerson’s  criterion  works  well  when  applied  to  problems  where  the  covariance 
matrix  is  such  that  the  above  problems  do  not  exist.  When  these  problems  arise,  as  in  Fig.  A1  or  in 
the  case  of  extraneous  interference  as  discussed  in  Section  5,  an  alternative  criterion  should  be 
applied. 


Fig.  A1  —  Interference  matrix  (17),  N  =  5 
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NORMALIZED  FREQUENCY 


Fig.  A2  —  Improvement  factors 
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Appendix  B 

PROOF  OF  ROBUSTNESS  THEOREM 


Define  the  vectors  S (0),  S,,S/by 


where  i  =  l,  2,  N  and  6  €  [0,2tt],  It  can  be  shown  that  the  sets  (Slt  S2 . SN]  and 

IS/,  S2  . . .  ,  S/yl  are  orthonormal  bases  for  complex  A/- space  CN.  This  additional  set  is  introduced 
for  convenience  in  the  proof  of  the  theorem. 

Lemma  1  —  Let  tv  €  CN  such  that  wTw*  =  /.  Then  1  <  max#  |  wrS(0)  |  2  <  N.  Equality  on 
the  right  is  achieved  if  w  a  S(-0max),  where  0^  is  any  0,  where  max#  |  wTS(6)  |  2  is  achieved. 
Equality  on  the  left  is  achieved  if  wTSx  =  wT S?  =  . . .  =  wT SN  =  l/'Z/V. 

Proof:  The  inequality  max#  |  wTS(6)  \ 2  <  N  is  argued  first.  Since  S, . SN  is  a  basis  for  CN, 

we  can  write 

*  =  E I  a<  1 2  =  1 

m  =  E  bj(6)Sh  £  \bi(9)\ 2  =  N.  (Bl) 

Now 


2 

E 
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but  by  the  Cauchy-Schwartz  inequality, 

r  "a  r 

=  £kl2  El*,!2 

L 1  J  L  j 


=  N 

with  equality  for  a  given  6  if  and  only  if  co  or  S(-0).  Since  |  wTS(0)  \ 2  is  a  continuous  function  over 
a  closed  interval,  the  maximum  exists  and  is  realized  at  some  0max  €  [0, 27r] .  Consequently  the 
preceding  inequality  implies  max^  |  h>  tS(6)  \2  <  N,  with  equality  if  w  a 

Next  it  is  proven  that  1  <  max^  |  wT S(6)  | 1 .  Since  (S/,  Sj  .  .  .  ,  S^j  is  an  orthonormal  basis, 
w  =  £  a, 5/  and  therefore  a,  =  h  ^S,  for  /  €  (1,2 . N\.  Observe  that  S(0,)  =  S,  for 

dj  =  27 ri/N,  from  which  it  follows  that  |  wrS(0,)  \  =  'fNai  and 

max  |  w  TS (6)  |  >  '^N  \  a ,  | 
o 

for  all  i.  Hence 

max|wrS(0)|  >  max  \^N  |a,  |  ,^V  |o,,  .  .  .  ,V/V  |  aN  \  ). 

6 

To  finish  the  argument  consider  two  cases:  (a)  |  a,  |  >  \/^N  for  all  i  6  (1.  2 . AM  and 

(b)  |  a,  |  <  1 J^N  for  some  /  €  ( 1 ,  2,  . .  .  ,  Ar).  In  case  (a)  the  result  follows  trivially.  In  case  (b) 

|  a,  |  >  1  nN  for  at  least  one  of  the  subscripts  /  €  (1,  2 . N] ;  otherwise  \  a,  \ 2  <  1, 

which  is  a  contradiction.  Therefore  the  maximum  on  the  right-hand  side  is  at  least  as  large  as  1. 
Hence 

max  |  w  TS(d)  (  >  1 

» 

and  consequently 

max  |  wTS(0)  \ 2  >  1. 
e 

Finally  equality  on  the  left  is  argued  if  a  j  =  a2  =  •••  =  fly  =  \/^N- 

max  |  wrS(0)  | 2  =  max  a,(S,')r  |  S(0) 

e  e  J 

r  r 

=  max  £(S/)r  £  bA^Sj 

V  1  J  V  J 
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where  b}  are  as  in  (a).  By  orthonormality  this  becomes. 


—  max 
N  8 


E  *;<*) 


Applying  the  Cauchy-Schwartz  inequality, 

<  -jj  max  E  I  bAe )  1 2  =  1  • 

The  result  follows  since  max#  |  wrS(0)  | 2  >  1. 

Proof  of  Robustness  Theorem: 

Using  the  definitions  of  Theorem  1  and  without  loss  of  generality,  we  assume  the  Filter  vector  w 
satisfies  wTw*  =  1.  Then, 

A(0)  =  10  log  Gc(6)  -  10  log  Gcn(d) 


=  10  log 


1  lwrS(fl))2 
y,  H' t Rw 


-  10  log 


1  1  wrS  (0)  1 2 

Jen  W  t/?  'W 


=  10  log 


Jen  Vf  t R  'h1 
y(.  H’ t Rw 


-  10  log  — 

lc 


1  + 


Q\ _ 1_ 

C  w  iRw 


(B2) 


S7s*  1 

where  the  input  SIRS  corresponding  to  the  covariance  matrices  R  and  R'  are  yc  =  - y  =  —  and 

No  i  <r 

S^S*  1 

ycn  —  - - - y  =  - - - — ,  respectively.  Note  that  G  is  a  continuous  function 


No ? 


1  + 


1 


1  + 


1 


-J  ~ 

N0  J  L  Alo¬ 
of  6  over  [0, 2-jt] .  Thus  the  maximum  of  G  is  obtained  for  some  6  in  the  interval  denoted  by  6^,  and 


^ max  max  Gc(0) 
e 


—  C?(0max) 


1  I  w  )  | 

yc  H'  iRw 
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Appendix  C 

COMPUTATION  OF  MMLE  FILTERS 


The  computation  of 


„t  S TR-'  S* 

min  C  R C  max - - - — 

C  [  «€A,  |  Cr  S  I  2 

and  the  corresponding  C  is  important  since,  to  the  authors  knowledge,  a  general  solution  has  not  been 
found.  However,  a  number  of  numerical  methods  are  available  to  compute  this.  Two  algorithms  are 
briefly  mentioned  here,  the  first  based  on  simulated  annealing,  the  second  on  a  deterministic  method 
called  the  directed  reduction  method. 


The  simulated  annealing  algorithm  [Cl]  involves  an  energy  function 

StR~'S* 


£|/|  =  min 
C| /I 


Cf/(+/?C|/l  max 


I11UA  T  , 

wa,  |  C|/|rS  | 2 


and  temperature. 


71/| 


7101 

log(l  +  i)  ' 


i  =  1.2,  . 


where  E\i  t  is  the  energy,  C(i)  is  the  filter  coefficient  vector,  and  T\i\  is  the  temperature  at  time  i. 
The  inner  maximization  over  0  can  be  approximated  by  uniformly  partitioning  0  6  [0,  2  7r)  and  discre¬ 
tizing  0.  Denoting  this  discretized  set  as  =  (0t  02 . 0s  I  function  becomes 


E\  1  j  =  min 
C|/| 


C|/|+/?C|/|  max 

0fA, 


ST  R~'S* 

i  ci»rs  i 2 


(Cl) 


For  the  examples  shown  in  this  report  S  was  between  seven  and  ten.  This  value,  however,  may 
need  to  be  increased  if  the  covariance  matrix  is  complicated  or  N  is  large.  New  coefficients  C'  are 
computed  as. 


c'  =  ciii  +  vi/i 

where  V\i\  is  a  vector  of  values  at  time  i  where  each  component,  at  each  time,  is  chosen  randomly  by 
a  zero  mean  IID  Gaussian  distribution  with  variance  proportional  to  7’|/|.  The  change  in  energy  is 
computed. 


A E  =  E\i  +  l|  -  £|i|. 
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If  A E  <  0  the  new  coefficients  are  accepted,  C(i  +  1)  = 
are  accepted  with  probability  p , 


C' 

(C  ')T(C')*' 


If  A E  >  0  the  coefficients 


-ae 


i.e.,  a  random  number  r  is  generated  uniformly  in  the  interval  [0,  1]  and  if  r  <  p  the  coefficient  vec- 

C' 

tor  is  C{i  +  1)  =  - — - ,  otherwise  C(i  +  lj  =  C{/(,  E\i  +  1|  =  E[i).  This  is  continued 

(C')r(C')* 

until  the  temperature  is  reduced  to  the  point  where  the  energy  reaches  steady  state.  Simulated  anneal¬ 
ing  algorithms  such  as  this  are  fairly  effective  and  easily  programmed.  Most  filters  designed  in  this 
report  took  between  10  and  60  min  of  running  time  to  approximately  reach  steady  state  on  a  Sun  3/60 
computer.  The  running  time  is  very  sensitive  to  the  initial  conditions  so  it  is  best  to  make  a  good 
guess  before  beginning.  Typically  filters  of  the  class  R~'S*  with  the  particular  0  €  h.A  that  minimize 
Eq.  (Cl)  are  useful  for  this  purpose. 


The  faster  directed  reduction  method  was  devised  as  described.  For  each  of  the  0,  € 
optimal  R~l S*  filter  coefficients  matched  at  that  frequency  are  predetermined.  Denote  these  coeffi¬ 
cients  as  fC|,  C2,  .  . .  ,  C§)  that  are  optimal  and  correspond  to  (0|,  02,  ....  0s)  respectively.  At 
each  step  i,  the  energy  is  computed  in  Eq.  (Cl).  Suppose  the  inner  maximum  occurs  at 
0  =  0/,/  €  (1,2,  ...  ,  Sj,  the  new  coefficients  are  then  determined  as 


C'  =  C{/)  +  kC, 


Cf /  +  1)  - 


C' 

(cY(cy' 


where  C /  is  the  optimal  weighting  corresponding  to  0/.  Note  that  the  overall  energy  Eq.  (Cl)  may 
increase  or  decrease  because  of  this  change.  It  should  be  monitored  to  see  if  steady  state  is  reached. 
This  algorithm  is  not  particularly  sensitive  to  the  initial  conditions  relating  to  the  time  to  reach  steady 
state,  assuming  the  algorithm  converges.  It  is,  on  the  other  hand,  less  stable  than  the  simulated 
annealing  algorithm,  and  the  initial  conditions  are  important  for  convergence  to  occur.  The  parameter 
k  controls  the  speed  of  convergence  and  is  related  to  the  algorithm's  stability.  If  the  value  is  too 
large,  the  algorithm  may  be  unstable.  For  the  applications  in  this  report  0.01  >  k  >  0.0001.  The 
filters  designed  for  this  report  took  between  1  and  15  min  of  running  time  to  approximately  reach 
steady  state  on  a  Sun  3/60  computer. 
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